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Abstract 

We propose, analyze, and implement interpolatory approximations and Filon- 
type cubature for efficient and accurate evaluation of a class of wideband generalized 
Fourier integrals on the sphere. The analysis includes derivation of (i) optimal order 
Sobolev norm error estimates for an explicit discrete Fourier transform type inter- 
polatory approximation of spherical functions; and (ii) a wavenumber explicit error 
estimate of the order 0{k~^N~^^), for i = 0,1, 2, where k is the wavenumber, N is 
the number of interpolation/cubature points on the sphere and depends on the 
smoothness of the integrand. Consequently, the cubature is robust for wideband 
(from very low to very high) frequencies and very efficient for highly-oscillatory in- 
tegrals because the quality of the high-order approximation (with respect to quadra- 
ture points) is further improved as the wavenumber increases. This property is a 
marked advantage compared to standard cubature that require at least ten points 
per wavelength per dimension and methods for which asymptotic convergence is 
known only with respect to the wavenumber subject to stable of computation of 
quadrature weights. Numerical results in this article demonstrate the optimal order 
accuracy of the interpolatory approximations and the wideband cubature. 
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1 Introduction 



The generalized Fourier integral operator J^, for a fixed phase function g and wavenumber 
K, G (0, oo), has the representation [H EI] 



Such integrals occur in various applications [H El [131 [IH ESI EDI EI] including simulation 
of waves scattered by the surface of an obstacle when an incident wave (say, a plane wave 
exp[iKX ■ d)) strikes the obstacle from direction d. In wave propagation models, for x E Q, 
the phase g{x) depends on the incident direction d and the observed point y E Q. Several 
fast algorithms and applications [H El |T3l [29l require accurate and efficient evaluation 
of such integrals for large number of observed points and also for wideband frequencies 
(that is, for n in (11. ip taking a wide range of values from very small to very large). 

Standard discretization techniques to evaluate such large numbers of generalized Fourier 
integrals is prohibitive, in particular for high-frequency wavenumber k and for Q C M^. 
This is mainly because standard quadrature rules require at least ten points per wave- 
length per dimension, even for low-order accuracy. This article is motivated by the recent 
work of Domfnguez et al. [9J on evaluating (11. ip in one dimension (with Q = [—1, 1] 
and g{s) = s on Q) and of Ganesh and Hawkins \T3\ on a fast algorithm to simulate 
high-frequency exterior acoustic scattering in three dimensions. 

In addition to the classical asymptotic methods [H EI], there is a large literature on 
efficient numerical evaluation of (II. ip . see [171123] and references therein. Such quadrature 
rules can be classified as Filon-type, Levin-type, or numerical steepest descent methods. A 
recent effort to avoid stability issues in such methods, for a one dimensional domain Q, is 
the shifted GMRES method [22] that requires solution of first order differential equations 
to evaluate I^f. The main aim in this literature is to prove asymptotic convergence 
with respect to the wavenumber k and address stability issues associated with computing 
weights for large k. We refer to the recent work [22^ i23j for details of the need to develop 
efficient methods to compute (II. ip even with Q = [—1,1] and discussions regarding the 
type of approach one needs to adopt depending on the properties of the phase function 
g. In particular, Filon-type methods are efficient if the phase function g does not have 
stationary points. 

For wave propagation applications, whether the phase function g has stationary points 
or not depends on the observation points. For solving the full obstacle scattering wave 
propagation models, observation points should be considered in all regions of the scatterer. 
Consequently, various types of techniques are required to evaluate (II. ip depending, for 
example, on whether the observation point is near or away from the shadow boundary and 
whether the corresponding phase function has stationary [H EI] or steepness [I3] points 
and also whether the density function / has singular points in Q. 

For a boundary integral equation reformulation of the wave propagation model exterior 
to a convex scatterer, efficient approaches to evaluate for stationary, steepness, and 
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singular point cases are discussed in fl8\ [13] (and references therein), respectively, for 
two and three dimensional cases. In such applications, Q in (11. ip is the boundary of 
the convex obstacle, which can be diffeomorphically mapped onto the unit sphere 

in M", n = 2,3. Such a global mapping property was used in [TBI [13] to transplant 
the wave propagation model to a boundary integral equation on §"'~^, n = 2,3 and 
linear combinations of appropriate polynomial basis functions are used to approximate 
the unknown surface current. 

The focus of the recent work of Domfnguez et al. [9] is to evaluate for a two dimensional 
scattering model [18] for the case when g does not have stationary points and / does 
not have singular points in Q. This leads to the requirement of evaluating (II. ip with 
Q = [—1,1] and g{s) = s, s E [—1,1]. We note that the integral considered in [9] is 
equivalent to (II. ip with = and g{x) = x ■ d, a; G when the incident wave 
direction is assumed, without loss of generality (due to rotational symmetry of §^), to be 
d = [0, 1]^. The main aim of this article is to carry out the three dimensional counterpart 
of the method and analysis in f9], by taking f2 = §^ and g{x) = x ■ d, a; G §^ in (II. ip 
with d= [0,0,1]^. 

For the sphere case, even basic tools that are used to accomplish an efficient method 
and analysis in p] is missing in the literature. A major contribution in this article is to 
first develop such tools and analysis, and hence derive an efficient Filon-type cubature 
and associated error estimates to evaluate the class of wideband integrals on the sphere 
defined by 

IJ= / f{x)exp{iKX ■ d) dx, d= [0,0,1]^, /«g(0,oo). (1.2) 

The Filon-type quadrature rules (and wavenumber explicit analysis) to evaluate the gen- 
eralized Fourier integral is based on the idea of replacing the density function / by an 
interpolatory approximation QnI G that can be efficiently computed. In addition for 
wavenumber explicit analysis, associated interpolatory approximation errors in Sobolev 
norms are required [9]. The choice of the finite dimensional space X^r, associated basis 
functions, interpolation points and simple representation of QnJ are crucial for efficient 
and stable approximations of the generalized Fourier integrals and analysis. 

For example, in [3], with X^- being the space of all polynomials of degree at most on 
[-1,1], 

4/ = y f{x) exp(i/€x) dx ^ I^^^f ■= j {QNf){x) exp(ira) dx, (1.3) 

where QnI interpolates / at + 1 (quadrature) points in [—1,1] that are Chebyshev 
points tj^N = cos(j7r/A^), j = 0, . . . , A^. In addition, QnI can be explicitly expressed using 
the discrete Fast Fourier Transform (DFFT) type linear combination of the Chebyshev 
polynomial basis functions T„, n = 0, . . . , N: 

N N 

QnI = 5^"(/,r„);vT„, (/,r„)^ = _^"cos(jn7r/A^)/(^.,7v), (1-4) 

n=0 j=0 
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where means the first and last terms in the sum are to be halved. In fll.4p . (-, ■)n is 
the discrete inner product quadrature approximation to the weighted inner product 

/■i dx r 

{f.9)--= f{x)g{x)-=== /(cos%(cos^)d^ 
J-i V 1 - Jo 

with a crucial property that 

(^l,^2)Ar = V^i, V'2 e Xat. (1.5) 

The interpolation operator with explicit DFFT-type representation and quadrature prop- 
erty in fll.4p - fll.5p does not require solving any matrix equation. Such an operator is known 
as the DFFT-type matrix- free interpolation operator [15]. The error analysis in [9j cru- 
cially depends on the classical error estimates of approximating 27r-periodic functions in 
the Sobolev space H^'-{-tt, tt), for = 0, 1: With ijjc{0) = ■?/^(cos6'), for < /x < i^, z/ > 1/2, 
there is a constant C^^^ such that [H] 

ll/c - iQNf)c\\H.(~.,.) < a,M^^^"1|/c|k^(-.,.), = 0, 1. (1.6) 

For the sphere case, with being the (A^ +1)^ dimensional space of all spherical poly- 
nomials of degree at most > 3, it is impossible to construct a DFFT-type matrix-free 
interpolation operator [25j. This seminal work of Sloan resulted in addressing several 
fundamental questions related to polynomial interpolation and numerical integration on 
the sphere, see [21 [16], [26] and extensive references therein. 

Using a new class of finite dimensional spaces Xjsf of spherical functions (some of which 
are not polynomials), Ganesh and Mhaskar [15] constructed matrix- free interpolation 
operators on the sphere. However, the analysis in [15] does not contain the Sobolev norm 
error estimate of the DFFT-type interpolation operators. In particular, for interpolation 
operator considered in this article we are not aware of error analysis in any norm. Deriving 
the Sobolev error estimates similar to fll.6p is one of the main contributions of this article. 

This paper is organized as follows. In the next section using a finite dimensional space in- 
troduced in [121 [U] we prove that the interpolation problem onto the space is well defined 
for any set of -|- 1 points in the latitudinal (elevation) angles in [0, vr] (that include the 
poles) provided that 2A^ equally spaced points are chosen in the longitudinal (azimuth) 
angle interval (— vr, tt]. In Section[3]we study Sobolev-like Fourier spaces that orthogonally 
decompose the standard Sobolev spaces on the sphere and establish properties of functions 
in such spaces. (Proofs of such technical results are deferred to Appendix [X]) In Section [H 
we prove optimal order convergence, similar to f ll.6p . of a DFFT-type matrix-free interpo- 
lation operator induced by the Gauss-Lobatto latitudinal quadrature points. (Major part 
of proofs in Section [Hare deferred to Appendix [B]) In Section [5] we introduce a Filon-type 
cubature to evaluate the wideband integrals fll.2p and prove that the associated error is of 
the order 0{k~^N~'^'), for £ = 0, 1, 2 where depends on the smoothness of the density 
function /. Numerical results in Section [6] demonstrate the robust theoretical analysis 
and efficiency of evaluating the wideband integrals for several low to high frequencies. 
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2 Interpolation on the sphere 



In this section we introduce a class of interpolatory approximations to continuous functions 
on Following the standard convention, we identify functions on the sphere in terms of 
functions defined using the spherical polar coordinates, 

F(^,0):= (F°op)(^,0), (2.1) 

where we used the natural spherical polar coordinates parameterization 

^ = 0) = (sin6'cos0, sin6'sin0, cos6')^, xeS"^. (2.2) 

Throughout this article, any function defined on the unit sphere §^ will be superscripted 
with "o". 

Using (12. 2p . we obtain the following periodic and refiectional symmetric properties of F 
in (EH) 

F(^,0) = F(0 + 27r,0) = F(^,0 + 27r), F(-^, + vr) = F(^, 0). (2.3) 

Let C(§^) be the space of all continuous functions on and define 

C := {F : ^ C I F satisfies (El]) for some F° G C(§2)}. (2.4) 
Observe that if F 6 C then 

F(0, ■), F{tt, ■) are constants. (2.5) 

Conversely, if F : — C is continuous and satisfies (12. 3 p and (12. 5 p then there exists a 
unique F° E C(§2) such that F = F° o p. Hence, C = C(§2). 

We first consider a finite dimensional subspace of C, depending on a parameter N, in 
which we shall seek interpolatory approximations. Such a space has been introduced and 
studied in [121 US]- We refer to details in [121 [IS] fo^' arriving at the space by restricting 
the space of all bivariate trigonometric polynomials of degree at N to those functions that 
satisfy the essential spherical function properties (12. 3p and (12. 5p . 

For M G N, we consider the finite dimensional spaces 

:= span! cos n6' : = 0, . . . , M} = spanj cos" 6* : ra = 0, ...,Mj, 

:= span! sin 61 cos ra6' : n = 0,...,M} (2.6) 

= span I sin ^ cos" 6' : n = 0,...,Mj. 

Throughout this article N eN with > 2. For p G there exist q G D^_2 and r G 
such that p{9) = r{6) + sin^ 6 q{6). Hence if p G D^^, we obtain 

p(0) = p(7r) = p{9) = sin^ 6 q{9), with q G D^_2. 
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Following [121 US]) we consider the finite dimensional subspace Xn C C: 

■= \poiO)+ ^ sin^ 6 pm{0) exp{im(j)) + ^ pm{0) exp{im(f)) 

^ -N<m<N -]V<m<JV 

even m^^O odd m 

Po e D^^, e e ^ 0, p2^+i G ©^_2} (2.7) 

= S ^ Pm(6')exp(im0) + ^ ^^(61) exp(im0) 

V ~N<m<N -N<m<N 
even m odd m 

P2e e ro^^r, P2£(0) = P2e{7i) = for £ 0, p2£+i e ©^_2}- (2.8) 
It is easy to see that the dimension of Af^r is 2N'^ — 2N + 2. 

Let = ^o<^i<''"<^Af = 7r with arbitrarily chosen set of — 1 latitudinal angles 
in (0,7r). In this article, we choose 2N equally spaced azimuthal angles (pk = kir/N, k = 
—N + 1, . . . ,N. These sets induce a discrete set of 2N'^ — 2N + 2 points on the unit sphere, 
by noticing from fl2.2p that p(^o, <Pk) and p{^n, (pk) are respectively the north and south 
pole, for = -A^ + 1, . . . , A^. 

Using these discrete sets of coordinates, we introduce the interpolation problem: For any 
F eC, find QnF e such that 

{QNF){^j,!f) = F{^^,^), j=0,...,N, k = -N + l,...,N. (2.9) 

For any F ^ C, since F(0, ■ ), F{n, ■ ) are constants, ( 12. 9p is equivalent to only 2N'^—2N+2 
interpolation conditions. Next we constructively show that the interpolation with the 
arbitrary — 1 elevation angles in (0, vr) is well-posed. Such a construction will also play 
a key role in developing an efficient approximation of the wideband integrals (11.21) . 

Proposition 2.1 For any F E C, the interpolation problem (\2.9\i has a unique solution. 

Proof. It is sufficient to show that such an interpolant can be always constructed. Set 

Fl:=F{^„!f), j=0,...,N, k = -N + 1, . . . , N, 
and define for j = 0, . . . , A^ 

-■= ^ H F^xp(-^), m = -Ar+l,...,Ar. (2.10) 

~N<k<N 

The vector P := {f^)m=-N+i,...,N C C^^ is indeed (a symmetric variant of) the result of 
applying the inverse finite Fourier transform to the vector F^^ := (F^)fc=_Ar+i_..._Ar C C^^. 
Hence, F^^ can be recovered from P via 

H= E /4exp(^), k = -N+l,...,N. (2.11) 

-N<m<N 
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Using (12. 3p . we derive 

F° = F(0,^) =F(0,0), = F{n,^)=F{7r,0), k = -N + I, . . . , N, 

and consequently 

/°=/^ = 0, forallm^O, f^ = F{0,-), /o^ = F(7r,-). (2.12) 

The interpolation problem: Find p2£ G such that 

P2i{Q = fL J = 0,...,iV, 

is equivalent, with the change of variables x = cos 6*, to a uniquely solvable polynomial 
interpolation problem on [—1, 1]. Moreover, fl2.12p implies that for i 0, 

P2t{0) = P2i{7l) = 0. (2.13) 

Similarly, for the odd coefficients, we solve the well-posed interpolation problem: Find 
P2i+i e such that 

P2i+i{Q = ^, j = i,...,iv-i, 

sm^j 

and set P2e+i '■= sin( • )p2e+i ^ '^°n-2- Hence we obtain 

P2e+i G 0^_2, P2e+i{Q = fii+i, j = 0, . . . , iV. 

If we define 

FNi6,(f)):= ^ p.mi9)exp{im(j)) e Xn, 

-N<m<N 

using (12. lip , we obtain 
M^j,w) = E Pm{QexpC-^)= J2 /4exp(^) = F^=F(0,t). 

-N<m<N -N<m<N 

That is, QnF = F^ is the unique solution of the interpolation problem (12. 9p . □ 

Note that for the well-posed interpolation problem (12.90 . a priori, we are free to choose 
the A^ — 1 nodes at the elevation angle. The quality of the interpolant depends crucially on 
the choice of the latitudinal interpolation points = < < ■ ■ ■ < = n . In particu- 
lar, for developing efficient Filon-type cubature to evaluate the wideband integrals (11.21) . 
we require an interpolation operator Qat on the sphere with optimal order convergence 
properties, similar to (II. 6p . in the Sobolev W norm for s = 0, 1. (Definitions of these 
norms are presented in next section.) 

That is, we require the Lebesgue constant of the interpolation operator Qa? to be uniformly 
bounded in T-L^ norm and grow only by 0{N) in the Ti^ norm. In addition for efficient 
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evaluation, similar to (11 ■4p . it is ideal to construct QnF with a DFFT-type explicit 
representation. 

The simple choice of ^j, j = 0, . . . ,N being equally spaced points in [0, vr] was considered 
in [111 [15] and the Lebesgue constant of the associated interpolation operator on the 
sphere in the uniform norm (that is, the standard norm in C) was proved in detail to be 
0{\og'^ N). The work in [15] in addition includes developing a grid-specific discrete inner 
product in A'jv and an associated orthonormal basis and hence a DFFT-type representation 
of Q^, see [151 Theorem 2.3]. While the space Xj^ in (12. 7p is same as that in [12] (with 
index m satisfying — < m < N) the finite dimensional subspace of C in [15] is slightly 
different (with index m satisfying \m\ < N). Results in [TJ] hold with |m| < replaced 
with -N <m< N. 

A DFFT-type matrix-free interpolation operator Qfj on the sphere with the + 1 lat- 
itudinal points based on the Gauss-Lobatto points was also introduced in [151 Theorem 
2.2]: Let {r^jjjLo be the nodes of the Gauss-Lobatto quadrature rule in [—1, 1]. That is, 
?7o = — l,?7Ar = l and ?7j for j = 1, . . . , — 1 are the roots of P^, where P/v is the Legendre 
polynomial of degree A^ |2j . The A^ + 1 latitudinal Gauss-Lobatto points are 

^. = Oj := arccos(r/^_j), j = 0, . . . , A^. (2.14) 

Thus Qff : C — )■ Xj^ is defined, for F G C, as the solution of 

{QiF)iOj,^) = Fie„^), j = 0,...,N, k = -N + l,...,N. (2.15) 

For a DFFT-type matrix-free analytical representation of Q^, we refer to [T5l Equa- 
tion (2.23)]. In this article we use Qfj for developing an efficient cubature rule for the 
wideband integrals. 

Unlike the case, convergence properties of Q|J- have not been established in any 
norm [15]. One of the key contributions in this article is to solve this open problem. In 
Section H] we will prove optimal order convergence properties of in the Sobolev "H* 
norm for s = 0, 1, after establishing various results in the next section. The convergence 
properties will be used in Section [5] to prove robust error estimates for a new class of 
cubature rules on the sphere, with respect to the wavenumber and quadrature points. 



3 Sobolev and Sobolev-like Fourier coefficient spaces 

The standard Sobolev spaces on the sphere 'H*(S^), can be introduced in several, but 
equivalent, ways. For instance, we can start from a smooth atlas for S^, i.e. local charts 
with a subordinated partition of unity, and define the space in terms of the appropriate 
Sobolev space of functionals on M?. The space so defined is independent of the chosen 
atlas used to describe the surface [H [2T] . 
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An alternative approach is based on using the eigensystem of the Laplace-Beltrami oper- 
ator on the sphere and construct a Hilbert scale which gives rise to the Sobolev spaces for 
any range. We follow the latter approach. For a sufficiently smooth function F° : E>^ ^ C, 
we consider the parameterized surface gradient and Laplace-Beltrami operators 



sin 6' a0 

1 d^F 
sin^ 6 dd)^ 



+ - 



1 d 



sin 9 89 . 



dF 

sm9—{9,<^) 



Corresponding to £^(§^), the space of all square integrable functions on the sphere, we 
define the equivalent space 

H°=|F:M2^C F satisfies (123])- (123]), ||F||^o < oo}, 

|F(^,0)|"sin^d^d0. (3.1) 



with 



n27r 



We have just used the fact that sin6'd6'd0 is the surface area element on the unit sphere. 
Similarly we define 

:= {F -.m^ \ F e H^, Vs2F G X H^] 

equipped with the norm 
1 



4 
+ 



|F||^o + ||Vs2F"2 



^0 



^0 



\F{9,(l))\^sin9d9 

2 



dF 
90' 



1 



sin 9 



d9 d0 + 





OF 


lo Jo 





sin^d^ d(f). (3.2) 



In order to introduce the equivalent high order Sobolev spaces through the eigensystem 
of the Laplace Beltrami operator, we consider the associated Legendre functions [21 121] 



(-ir 

2"n! 



^m+n 



dx'^+" 



(3.3) 



Define 



, ) ^ : /^-l(cos^), Q„-"^:=Q^ 0<m<n, n = 0, 1, . . . 
2 [n + my. J 

(3.4) 

Gm{4>) '■= I exp(im0), m E Z. (3.5) 



and set 



The spherical harmonics [21 [21], polynomials of degree n on the sphere, can be defined as 
y7(^,0):=(-l)('"+H)/2g;n(^)^^(^)^ m = -n,...,n, n = 0, 1, . . . (3.6) 

It is well known that {Y^ : n = 0, 1, 2, . . . , |m| < n} is an orthonormal basis of 
^ £2(g2) g^^^ [2X1 (2.4.175)] 

/))sin^d^d0 = 



^0 



n(n + 1), if {n, m) = {n\ m'), 
0, otherwise. 



(3.7) 

The above result follows from the fact that {Y^} are eigenvectors of A§2 [211 (2.4.23)]: 

-As2Y;^={n{n + l))Y^, m e Z, n = \m\, . . . , oo. (3.8) 
Thus {(n + l/2)^^y^ I 77, = 0, 1, 2, ... , |m| < n} is an orthonormal basis for Ti^. 

For any s G M we introduce the norm 



oo n 



n=0 m=— n 



where 



7r /'27r 



^0 



F(^,0)r™(e,0) sin^d^ d(/). 



The Sobolev spaces "H** for s > are then defined as 

H' := {F G H° I \\F\\hs < oo}. 



(3.9) 



(3.10) 



Observe that for any F 



oo n 



m 
n 1 



n=Q m=—n 



with the series converging now in "H*. Thus, spanjy^ : n = 0, 1,2, . . . 
dense in 'H'* for any s. 



(3.11) 

|m| < n} is 



with We point out that for s = 0, 1, the norms in (13. 9 p and in fl3.ip - fl3.2p coincide. 

Moreover, 



|2 



\F\inj2 



|A§2F||^o + 7l|V§2F||^oxwo + — ]|F]|^o 



TV p27r 

Jo 



1 d^F 



16' 

1 d f . JF^ 
sm 9 o9 \ o9 



+ 



4 



9^' 



1 



sin ^ 



d^ d(p + 



'0 ^0 

+ \Fi9,<P)\hm9d9d^. 



orthonormal with respect to this inner product. 



1 /"Stt 


OF 


4/0 I 






= 0,1,2, 



sin 9 d9 d0 

2 

sin^d^ d</) 

(3.12) 
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3.1 Sobolev-like Fourier coefficient spaces 



For s > 0, we introduce a decomposition of W as an orthogonal direct sum. Such a 
decomposition, roughly speaking, comes from interchanging the order in which the series 
in (13. lip is summed. In order to simplify the exposition, we start first with F G H". 
Hence 



oo n oo X oo \ 

^(^''^) = E E KmY:riOA)= E ( E KMe)^^ 

n=0 m=—n m=— oo ^-o— i™i / 

oo 



■ n=|m| 



E [J'mFme^ict)). 



(3.13) 



The function F{9, ■ ) is almost everywhere a function in L^(0,27r) and as such it can be 
expanded in its Fourier series. By matching the Fourier coefficients of the same order, we 
conclude that J^mF admits the integral expression 



T F 



27r 



F( -,0)6^(0) d0 



2-K 



F(-,0)e_„(0) d0. 



Furthermore, 



I^IIho 



2-K 



|F(^,0)|M0 



sm 



^d^^= E ii-^-^i 



2 

sin 



m=— OO 



with 



|/(^)|2sin^d^ 



L "'0 



1/2 



Then, if is the corresponding weighted space, we obtain 

oo 

= {f® I /l(0,.) e /(■ + 27r) = /,/(-■) = (-1)"^/}, 



m=— oo 



where 



(/®ej(0,</.) :=/(^)e. 



f{6) exp(im(; 



Observe that the sum above is actually orthogonal. 



(3.14) 



(3.15) 



(3.16) 



We extend this decomposition to l-C^ by introducing new spaces, one for each Fourier 
mode. Let 



' ' m 

endowed with the image norm 



(3.17) 



11/ ®e 



m II "H^ 



(3.18) 
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Clearly, = Wt^. We note that the definition of imposes naturally some period- 
icity and parity conditions on its elements, namely 

feW^ =^ /(.+27r) = /, /(-.) = (-ir/ (3.19) 

which are simply refiections on those satisfied by the elements of W (see fl2.3p ). 

It is easy to see that, for each m G Z, is a Hilbert space and a complete orthonormal 
set in is {(n + 1/2)-"Q;^ : n > \m\}. Hence, 



oo / oo \ 1/2 

n=|m| n=|m| 



with 



fm{n) := / f{e)Q^{e)sinede= / (/®e„)(^^,0)y;-(^,0)d^d0 
Jo Jo Jo 

= (/ ® e„)„_^. 

Clearly, : 'H'^ is continuous and onto. Moreover, using fl3.13p - fl3.15p . 

oo oo 
m=— oo m=— oo 

with the first series converging in W. That is, we obtain an orthogonal direct sum 
decomposition of 'H'^: 

oo 
m=—oo 



3.2 Properties of Fourier coefficient spaces 

We derive some properties of that we will use in this article. Proofs of the results in 
this section are given in Appendix |Al For further study of such spaces, we refer to flU\ 
where a closely related family of spaces, namely that obtained by making the change of 
variables x = cos6, are studied. 

It is easy to see that for r > s the injection C is compact. Moreover, because of 
(I320D, 

\\f\\w^<{\m\ + lr^\\f\\w^, \/r>s. (3.22) 
Equations f l3.2ip and fl3.22p prove that for r > s 

oo 

E (1^1 + iy^''~'^\\J''-n'w^ < Wn'nr. (3.23) 
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Hence, 



-N<m<N 



-N 



|2 



m=N+l m=— oo 



1/2 



< (iV + i)^-'-[ 5^ (|m| + i)2('-^)||j-„F 

|m|>Af 

< (iV+i)^-'-||F||^. 



for all r > s. 



i2 



1/2 



(3.24) 



Proposition 3.1 For all m G Z, 



2 



1 
4 



d0 
sin 



Moreover, if \m\ > 2, 



2 



\f{d)\^ sinddd. (3.25) 



+ irri^ - Arn^] 



sin^ ^ 



+ ^ r|/'(^)psin^de+(l + 2m2) r\f'(9)\'^+ F \f"{9)\' sin 9 d9. 

(3.26) 



The norms f l3.25l) - p.26l) contain three and six integral terms respectively. We can, how- 
ever, define equivalent norms with fewer terms. Let 



2 



2 



r\f{e)\'^+ r\f'{e)\\m9d9 
Jo sm 9 Jq 

|2 2 I ////)M2 , / I W///1M2 



(3.27) 







sin^6' 



+ I If 





. + / |r(e)|^sin0de. (3.28) 
sm Jo 



The following result shows that, for i = 1,2, the spaces and are equivalent for all 
|m| > i with equivalence constants independent of m. 



Corollary 3.2 For all m E Tj with \m\ > 1, 



w,} < 



< ^ii/iu^ < 4 



(3.29) 



Moreover, for all m E Z with \m\ > 2, 



1 

7! 



Z2 



< ^WfWwi- 



(3.30) 
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Remark 3.3 Notice that if / G = Zj^ for m ^ 0, then and its derivative are 
in L^(0,7r). Thus, from the Sobolev embedding theorem |/p is continuous, and so is /. 
Moreover, / vanishes at the end points {0, tt}. This regularity property need not be true 

1 1 /2 

for functions in W^. For instance, | log | sin^| | e Wq but this function is obviously not 
continuous. □ 



We complete this section by studying the regularity of from a classical Sobolev point 
of view. Let us introduce the 27r-periodic Sobolev spaces 



Hi-' 



loci 



/ = /(- + 27r)} 



(3.31) 



endowed with the norm 



^^:=l/(o)r + EH'n/Mr, Km)^^r f{e)eM--^me)de. (3.32) 



For r = 0, II • 11^0 is simply the L^(0, 27r) norm, and for non-negative integer values of r, 

# 

an equivalent norm is given by 



Jo 



1/2 



(3.33) 



Proposition 3.4 For all r > there exists > independent of m and f such that 



Further, 



(3.34) 



< C(l+|m|)||/||^i, 



(3.35) 
(3.36) 



with C independent of f and m. 



Proposition 3.5 There exists C > such that for all f e Wq 



r\f{e)\'d0 < c\r\f{9)\ 

Jo IJo 



'sin^d^ 



+ 1 1^ \ fie)]' sine de 



1/2 



|/'(^)|2sin^d^ 



1/2. 
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4 Optimal order convergence analysis 



In this section, we develop several results and hence prove the following optimal order 
convergence property of the interpolatory operator Qfj on the sphere. 

Theorem 4.1 Let s = 0, 1. For all r > s + 3/2 there exists Cr > such that for all 
F G W the following optimal order estimate holds: 

(1 - s)\\ Q^F - Fll^o + sN\\ QiF - < a,N-'-\\FUr. (4.1) 

We demonstrate the optimal order convergence property with several numerical experi- 
ments in Section [6l 

The discrete spaces and the functional frame where the analysis will be carried out have 
been already introduced in previous sections. The proof of Theorem 14 . II will rely on using 
certain semi-discrete operators and two interpolation operators and in the polar 
angle 9. These operators facilitate a representation of the interpolation error as a sum of 
that introduced in the approximation of the central 2N Fourier coefficients of F plus the 
error arising from ignoring the tail of the Fourier series (I3.2ip . 

4.1 Fourier analysis 

Recalling the definition of in (13.131) . the aim of this part is to investigate properties 
of the semi-discrete maps 

oo 

(p^F) {6) := (-^m+2£7vF) (0), m = -N + l,...,N. (4.2) 

l=-oo 

The following bound will be used repeatedly in this article: For r > 1 and for all |m| < 

oo 

e^o 1=1 
Lemma 4.2 For all r > 1 there exists Cr > such that for any F G 

oo oo ^ oo oo 

m=—oo n=\m\ ^ m=—oo n=\m\ 

Moreover, 

F(0, ■) = (-^oi^)(0), F(7r, ■)= (•^oi^)(vr), 

and therefore F & C. 
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Proof. Since F E "H^^^, for some e > 0, using fl3.6p . 



oo oo 



oo oo 



m=—oo n=\m\ m=—oo n=\m\ 



< 



v ^ m=— oo n=lml n=0 m=—n ' 



■ m=-oo n=|m| 
oo n 



2+2£ 



'- n=0 m=— n 



1/2 



2e 



^ n=0 



1/2 



<aiiFi 



In the penultimate step, we have used the addition theorem [2T| Theorem 2.4.5] 

n 

r?i=— n 

In particular, F is continuous. Moreover, for m 7^ 0, (5^(0) = <55^(vr) = for all n. Hence 
using fl3.13p . 

00 

F(0, 0) = 5^ ^n,oQ° (0) = ( J-qF) (0), V0. 

n=0 

Analogously, it can be shown that F(7r, ■ ) = [TqF) (tt). Thus F EC. □ 



Lemma 4.3 For all\m\ < N and r > 1, p|^-F G C[0, tt] and 

(p^F) (0) = (vr) = 0, ^/ m ^ 0, (4.4a) 

{p%F){0) = (j-oF)(O) = F(O,0), (pMW = {J'oF){n) = F{7i,4>), V0 G M. (4.4b) 



Proof. Follows from the definition of J-'m in f l3.13p and Lemma 14.2 



□ 



Given a general interpolation operator Qat : C — )■ Xn fulfilling the conditions stated 
in (12. 9p . we consider the associated one dimensional interpolants: For / G C[0,7r] find 
q^/ G and q%f E D^_2 such that 



^%fi^^) = /te), t = 0,...,N, 
c&fiQ = /(O), J = l,...,iV-L 



(4.5a) 
(4.5b) 



With the help of these interpolants we are able to write QnF in terms of the Fourier 
coefficients of F. 



Lemma 4.4 For all F E W with r > 1 

(QivF)(^,0)= (qM^)Wem(0)+ (qM^)Wem(0). (4.6) 



-iV+l<m<iV 
even m 



'N + l<m<N 
odd m 
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Proof. Note that, because of (14.41) . for even m 7^ 0, {(1%Pn){^) = (QArPjv)(^) = which 
proves that the element of the right hand side of (14. 6 p is contained in X^. We then just 
have to prove that the right hand side of (14.61) satisfies the interpolation conditions (12. 9p . 

Since F G T-T, with r > 1, Lemma 14.21 shows that the Laplace series of F converges 
absolutely. Hence the following manipulations are fully justified. For k = —N + 1, . . . ,N, 

00 00 

m=-oo ~N+l<m<N e=-oo 

= J2 (fl {J'm+2eNF)ie))emm 

-N+l<m<N ^£=-00 ^ 

E iPNFmcr^m. 

-N+l<m<N 



Therefore, for j = 0, . . . , N and k = —N + 1, . . . , N , 

nivm) = E iPNFm)e^{i§) 

~N+l<m<N 



-N + l<m<N -N + l<m<N 

even m odd m 



Thus we obtain the representation (14. 6p □ 



Denote by 

q%^^,:C[0,n]^B%, and q%^^, : C[0,n] ^ (4.7) 

respectively the interpolant projections defined in (14. 5 p with nodes {9j}^^Q, the Gauss- 
Lobatto points. Then using Lemma [4.41 

QiF{d,<p)= E {^NSPNF){e)eU<P)+ E {^NsPNF){e)eU<P). (4.8) 

-N + l<m<N -N + l<m<N 

even m odd m 

For any m G Z, we introduce a simpler notation 

._ / q^,gi, if rn is even, 

■ I q^,gi, if ^ is odd. ^ 

We recall the spaces Z^, for i = 1, 2, induced respectively by the norms (I3.27P and (13.281) 
and that for |m| > i, these are equivalent to the spaces in (13.170 with respective norms 
given by (I3.25P and (13.260 . The next two results give properties of q^gj when applied 
to elements in these spaces. Derivation of these results require new technical tools and 
hence we postpone the proofs of the following two results to Appendix [Bl 
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Lemma 4.5 There exists C > such that for all f G Zj^ with m ^ and N > 2 

l|q^,gi/IUL<^'""" ^ 
Besides, for all f G with m>2 

\\^nJ-\W^<c 

where C is independent of m and N . 
Proof. See Lemma [B .41 and Corollary IB. 91 



(4.10) 



(4.11) 



□ 



Proposition 4.6 Let s = 0, 1. Then, for s = and r > 1 and for s = 1 and r >2, there 
exists Cr- > such that for any N > 2 and f G WJ^ 

llq^,gi/-/lk^<aiv^"1l/lk^- 

Moreover, for m 7^ and s = the result holds also for r = 1. 



Proof. See Propositions IB.6 1 and IB.lOl 



□ 



The main result of this section is stated and proven below. 

Theorem 4.7 Let s = 0, 1. For r > 3/2 + s, there exists Cr > such that for all N >2 
and \m\ < N , 



I^^rXPnP ~ ^mF\\w^ < CrN 



m,+2lNF\\'^r 



m+2tN 



e=-oo 



1/2 



< aA^^-^||F||^..(4.12) 



Proof. We have 

sm sin 

For the first term in the bound, we obtain 



2 

sin 



l^ll 



< 



2 



< 2aiv-2r 



\2\m + 2m\^'-^\\J'^+2mF\\l2 + N-^J'm+2mF\\l^ ) 



^ |m + 2£iV|2^^-2(|m + 2iN\-^ + iV-^) || J-^+s^atFU'h/i _ J 



m + 2lN' 



m+2£NF\\'^r 



(4.13) 
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where we have apphed successively (14. 3 p , fl4.10p of Lemma I4.5[ Corollary 13.21 and fl3.22p . 
On the other hand, Proposition 14.61 yields 



The result for the case s = follows now from fl4.13p and fl4.14p . 
For s = 1, we proceed as before, but using now (14.1 ip of Lemma H75l 



(4.14) 



l-E'illv^i — 



^ |m + 2m\ 



4-2r 



£^0 



^|m + 2£iV|2'^-iq-J-^+2£^F| 



2 



Y,\m + 2lN\^'-W\J^m+2mF\ 



< CrN' 



■4-2r 



1^0 



2 

iy2 



m+2tN' 



< ICrN 



2-2r 



m+2lNF\\'w 



m+2lN 



(4.15) 



Finally, Proposition 14.61 yields 



|-E'2||wi < CrN^ ^WJ'rnFWwi, 



□ 



4.2 Proof of Theorem 14.11 

For s = 0, 1, using (14. 8 p and (I3.2ip we obtain 



(-N oo X 

-iy^7,t:^iv m=-oo m=Af+l^ 

+ D2. 



The last term can be estimated using (I3.24p 

D2 < N^'-^'^WFWnr. 

On the other hand, and with the help of Theorem 14. 7[ we deduce that 

00 

-N+l<m<N l=-oo 
00 

|2 



< aA^2s-2r J2 W^mFfwr^ = CrN^'-^'-\\F\ 



m=— 00 



□ 
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5 Cubature and analysis for wideband integrals 



In this section, using the tools developed in previous sections, we propose efficient ap- 
proximations to the class of wideband integrals on the sphere defined in fll.2p and analyze 
the error in such cubature approximations. 

In particular, for F G T-C with r > 1 and / = F° in f ll.2p . the wideband integrals (with 
wavenumber k G (0, oo) and incident direction d = [0, 0, 1]"^) can be written as 

I^{F) := / / F(0, 0)exp(iK cos sine d^d^ = / F°(a;) exp(ifi;a; • d) da;. (5.1) 
Jo Jo is2 

We are interested in robust cubature rules that give accurate results for small, large and 
very large values of the wavenumber k. 

Following the ideas of Filon rules, or those more general product integration rules (see for 
instance |I9]), we propose 

4(F) ~ / / ( Qf^F) (e,0)exp(i«; COS e) sine ded0=: /J^(F). (5.2) 
Jo Jo 

Using f l3.14p and Lemma [4. 4^ we obtain for F G with r > 1, 

4(F) = V2^ (JqF) (e)exp(iK cose) sine de, (5.3) 
Jo 

/g|^(F) = V2^ (Jo Q^F)(e)exp(iK cose) sine de 
Jo 

= V2n I (q^giP^F)(e) exp(ifi; cose) sine de. (5.4) 
Jo 



Consequently, 

l4(F)-/J^(F)|<v^ 



( J'oF - q^,giP^F) (e) exp(iK cos 6) sin 6 dO 



Before entering into the analysis of such approximations, we observe that computation 
of q^p^F does not require evaluations of either J^21nF for all £ or Q%F. This is a 
substantial computational advantage even compared to using the matrix-free FFT type 
representation of Q%F. An efficient algorithm to compute the cubature, based on the 
proof of Proposition 12.11 is as follows. 



Algorithm 

1. Set 

/° = F(0,0), /o^ = F(7r,0), fo:=^ E ^(0, f)' J = 1, • • • , AT - 1. 

-N+l<k<N 
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2. Solve the one dimensional interpolation problem at Gauss-Lobatto based nodes: 

]D)%3PN s.t. p^(e,) = /J, j=0,...,iV. 

(Observe that Pn = ^qAr,giPAr-F-) 

3. Compute 

2vr / j9iv(^) exp(iKcos^) sin^d6'. (5.5) 
Jo 



Now we require a robust approach for evaluating (I5.5p . After performing the change of 
variables 6 = arccos x we have to evaluate 

1 

PAr(arccosa;) exp(iKa;) dx. 

-1 

Since pAr(arccos ■) is now a polynomial, there are some accurate, robust, and fast methods 
in the literature for evaluating this integral, see [9] and references therein. In our numerical 
implementation, we use the approach developed in [H] for evaluating the one dimensional 
integrals fl5.5p . 

Crucial to our wavenumber explicit error analysis of the cubature rule is Theorem 14.71 
which is essentially the main result for proving optimal order error bounds for the inter- 
polatory approximation Qf^F in the Ti^ and norms. In addition, for the error analysis 
we require additional convergence estimates in stronger norms, defined fl3.3ip - fl3.33p . 

Lemma 5.1 Let g G and qn G D^. Then for all s > and r > s + 1/2 there exists 
Cr > so that 

\\g - ^Tvbj < Cr [N'-'^'/^g\\ws + N'^'/^g - g^W^o + N'-'/^g - g^Ww^] . (5.6) 



Proof. We consider the following projection on 



"AT- 



N 



SNg{0) '■= ''^^ g{m) cos m9, g{m) :- 



m=0 



J-j:^g{9)d9, ifm = 0, 

^ J'^^g{9)cosm6d9, if m > 0. 



Hence for even g, as it is in our case, it is easy to prove from f l3.32p (see also j24]) that 

bNg-gWni < N'llgWw^- 

The other result we need in this proof is the well-known and easy to prove inverse in- 
equality 

\\gN\\H^^ < N^^'MgNlln^, ^gN e span {exp(im ■ ) : < ^} ' s>r. 
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Using the above inequalities and fl3.34p . we obtain 



19 - 9n\\h^^ < Wg -SNg\\H^ + \\sNg - gNliH^ 



< iV^-'-+i/l(?||^._v2 + N'WsNg - gN\\H% 

< iV-^+i/2||^||^.-v2 + N^iWsNg - g\\Hi + \\g - gNWnoJ 

< CrN-'-^-^'/^WgWws + V2N'\\g - gNh^o,.)- 
We next apply Proposition 13.51 to the last term and obtain 

1 /2 

\\g - gN\\LHo,n) < C[\\g - gN\\l2, + \\g - gNh^. Ug - gNYh^. ] 

'- sm sin sin-i 

< C{l + N'/')\\g-gr,h. + N-'/^g - g^yW^. . (5.7) 
The result (15. 6p now follows by applying Proposition 13. II to the right hand side above. □ 



Now we are ready to prove the wavenumber explicit error bounds for the cubature ap- 
proximations lfj^{F) of the wideband integrals Ik{F). Such error bounds provide robust 
error estimates for convergence with respect to the number of cubature points and the 
wavenumber. 

Theorem 5.2 For i = 0,1 and r > 3/2 or i = 2 and r > 4 there exists Cr such that for 
all F , 

\UF) - lt,Nm < CrK~'N^^'^-^\\F\\nr, (5.8) 

where 



( 0, I 


fi = 


V{i):={ 3/2, , 


fi=l 


( 4, t 


fi = 2 



Proof. Let 



rN := V27r(q^^,giP^F- J-qF). 



For F smooth enough (precise information about the regularity requirements on F will 
be given below), r^(0) = 'r^(7r) = and therefore the function 



^ if^G(0,vr), 
sm 6 

r^(0), if^ = 0, 

-r^(7r), if^ = 7r. 
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is continuous. Hence, using integration by parts we derive 



IK 



IK 



r^^O) exp(ifi;cos6') sin 6* (16* 



r'^{6) exp(iK cos ^) 



(Pn{(^) exp(iKcos^) sin^d6' 



(5.9) 
(5.10) 



;'{)n{G) exp(iKcos6') 



1 



=0 K-^ 



V9'jv(6') exp(i/tcos6') d^^ . 

(5.11) 



Hence, (15. 9p implies that 

\In{F) - h,N{F)\ < ¥n\\l^. 

I \ / 1 \ y \ II II gjjj 

whereas (15.101) yields 

Finally, (15. lip shows that 

\h{F) - h,N{F)\ < llfNm + |¥^7v(vr)| + r |^'^(^)| d^' 

L Jo 

Proceeding as in the proof of Theorem 2.2 in [9], we obtain 



(5.12) 
(5.13) 



|<^;v(0)| + |<^(7r)| <C||r^||^,3 



and that for all £ > 0, 



|<^'^(^)| d^ < CeiN'/^WrNhno,.) + iV-^||rAr||^7/2+.)- 

# 

Hence, for all s' > 7/2 we deduce 



(5.14) 



Using fl5:T2|) and fl4:T2D of Theorem O with s = yields 

IkAflUl^ = lk7v||H/o < CrN'''\\F\\nr. 

Thus dEHD holds for £ = 0. 

For i = 1,2, we need to bound ||rAr||/^s for different values of s. Applying (15. 6p . with 
g = J-qF and gjy = qJ^p^F and (I4.12p of Theorem 14.71 we obtain for r > s + 1/2 



< Cr 



e=-oo 



1/2 



(5.15) 
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Applying flS.lSp first in (I5.13P with s = 1 and second in f l5.14p with s = 3, s = and 
s = s' < r - 1/2, we obtain the result (15. 8p for £=1,2. □ 

6 Numerical Experiments 



In the section, we demonstrate the algorithms and analysis developed in this article for in- 
terpolatory approximations of functions on the sphere with various smoothness properties 
and efficient evaluation of the wideband integrals in ( 15. ip induced by such functions. For a 
fixed observation point x* G and a smoothness parameter s, we consider F° : — ^ M, 
defined by 

F:{x) = \x-x*\'= [{xi - xlf + {X2 - xlY + (X3 - xlf] xe §^ (6.1) 

where x = [xi,X2,X3]^ and x* = [xl,X2,xl]'^ . 

Clearly, for s < 0, F° is a discontinuous spherical function; for s > and even, F° is 
a spherical polynomial; and for s > 1 and odd F° e ?^''(§^), for all r < s + 1. For 
numerical experiments, we chose x* = [2/3,1/3,2/3]'^; s = 1,3,5 and various wideband 
wavenumbers n = 10™ with m = —2, —1, 0, 1, 2, 3, 4, 5. 

For s = 1, 3, 5, Theorem 14.11 indicates that the estimated order of convergence (EOC) 
for approximating Fg by QfiFs in the T-L^ and Ti^ norms are approximately s + 1 and s. 
Results in Tables [1], [21 and [3] substantiate the efficient approximations and analysis for 
spherical functions with various smoothness properties. The Sobolev norm errors were 
computed by approximating integrals in (13. ip and (13.21) using rectangular rules with over 
160, 000 quadrature points. 

Further, for s = 1, 3, 5, using Theorem 15.21 the error in approximating the wideband inte- 



grals by the new class of cubature. 



is similar to that of || Q%Fs — Fs\\y^o 



UFs)-ltNiFs 

for all wideband wavenumbers and that the cubature error is much smaller for large 
wavenumbers. Results in Tables HI O and [6] demonstrate the very efficient approach de- 
veloped in the article to approximate the wideband integrals and associated error analysis. 
For a fixed and for large k > 100, results in Tables HI [5|, and |6] indicate O^k,^"^) decay 
of the cubature error. This is the case even for Fi that does not satisfy the smoothness 
assumption in Theorem 15. 2[ In particular, for a fixed large k, accurate solutions have 
been obtained even for small values and hence we do not observe consistent error im- 
provement for certain increment in values of A^ for the Fi case. Consequently, similar 
to [U Table 2] and as explained in [9], Experiement 2], for a fixed large k, the estimate 
in Theorem 15.21 for Fi is only an upper bound for the error. 

Using the algorithm described in Section 0, in addition to obtaining high-order accuracy, 
the total MATLAB computation time taken to produce all results in TablesHJEl andlHlwas 
less than one-tenth of one second of a single core of a 3.4GHz INTEL CORE 17-2600 CPU. 
This is a marked computational advantage over (i) standard cubature approximations that 
require at least ten points per wavelength to compute such wideband integrals with low- 
order accuracy; and (ii) asymptotic based approximations and analysis that are applicable 
only for very large wavenumbers. 
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Table 1: V,^ and T-L^ errors in approximation of Fi by Qf^Fi and EOC 



N 


\\Q%F^ - F^\\no 




\\Q%F^ ^ F,\\n-. 




5 


4.0077C-02 




3.3636E-01 




10 


1.1383e-02 


1.8159 


1.6166e-01 


1.0570 


20 


2.7244e-03 


2.0629 


8.9315e-02 


0.8560 


40 


6.7527e-04 


2.0124 


3.9452e-02 


1.1788 


80 


1.6554e-04 


2.0283 


2.0598e-02 


0.9376 



Table 2: Ti^ and "H^ errors in approximation of F^ by Qf^i^s and EOC 



N 








EOC{n') 


5 


9.1450e-03 




6.4589e-02 




10 


5.7752e-04 


3.9850 


6.8649e-03 


3.2340 


20 


3.4725e-05 


4.0558 


9.7395e-04 


2.8173 


40 


2.2636e-06 


3.9393 


1.0621e-04 


3.1970 


80 


1.4386e-07 


3.9758 


1.4652e-05 


2.8577 



Table 3: and Ti^ errors in approximation of F5 by Q^F^ and EOC 



N 


1 Qn^5 — FrjWno 


EOC{n°) 


1 Qn^^ — FrjWn^ 


EOC{H^) 


5 


8.7351e-03 




5.9412E-02 




10 


1.0756e-04 


6.3436 


1.2333e-03 


5.5901 


20 


1.5172e-06 


6.1476 


4.0722e-05 


4.9206 


40 


2.5432e-08 


5.8986 


1.1236e-06 


5.1796 


80 


4.0461e-10 


5.9740 


3.9233e-08 


4.8399 



Table 4: Cubature errors 



for K = 10^ 



m 



-2.-1.0. 1.2.3.4.5 



N 


K = 10-2 


K = 10-1 


K = 10° 


K = lOl 


K = 102 


K = 10^ 


K = 10^ 


K = 10^^ 


5 
10 
20 
40 


2.0e-03 
1.3e-03 
3.7e-05 
1.6e-05 


2.1e-03 
1.3e-03 
3.7e-05 
1.6e-05 


3.2e-03 
1.3e-03 
3.8e-05 
1.6e-05 


l.le-02 
3.6e-03 
1.3e-04 
1.7e-05 


7.3e-05 
2.6e-05 
1.4e-05 
1.4e-05 


9.5e-07 
1.8e-07 
1.4e-07 
8.1e-09 


7.7e-09 
1.4e-09 
1.8e-09 
9.7e-ll 


7.4e-ll 
1.3e-ll 
1.9e-ll 
9.7e-13 


Table 5: Cubature errors 


UFs) - itAFs) 


for K = 10™, m = -2, -1, 0, 1, 2, 3, 4, 5 


N 


K = 10-2 


K = 10-1 


K = 10° 


K = lOl 


K = 102 


K = 10^ 


K = 10^ 


K, = 10^ 


5 

10 

20 
40 


1.2e-04 
3.0e-05 
1.7e-07 
2.3e-08 


1.4e-04 
3.0e-05 
1.7e-07 
2.3e-08 


3.9e-04 
3.0e-05 
2.0e-07 
2.3e-08 


1.7e-03 
2.3e-04 
1.2e-06 
2.8e-08 


1.5e-05 
1.7e-06 
3.1e-07 
2.8e-08 


1.9e-07 
2.1e-08 
3.2e-09 
2.1e-ll 


1.6e-09 
1.7e-10 
4.0e-ll 
3.2e-13 


1.4e-ll 
1.6e-12 
4.1e-13 
3.3e-15 


Table 6: Cubature errors 


UF,) - itr^iF,) 


for K = 10", m = -2, -1, 0, 1, 2, 3, 4, 5 




K = 10-2 


K = 10-1 


K = 10° 


K = lOl 


K = 102 


K = 10^ 


K = 10^ 


K = 10^ 


5 

10 
20 
40 


3.1e-05 
1.7e-06 
1.4e-09 
9.0e-ll 


7.0e-05 
1.7e-06 
1.4e-09 
9.0e-ll 


1.5e-04 
1.8e-06 
2.4e-09 
9.0e-ll 


8.1e-04 
4.5e-05 
2.9e-08 
1.3e-10 


7.2e-06 
4.2e-07 
1.9e-08 
2.7e-10 


l.Oe-07 
5.5e-09 
1.9e-10 
8.1e-13 


7.0e-10 
4.3e-ll 
2.3e-12 
3.2e-15 


6.3e-12 
4.1e-13 
2.3e-14 
4.8e-18 
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A Proofs of results in Section 



K2 



Proof. (Proposition 13.11 ) To establish this result we will use fl3.18p with (13. 2p and 
(I3.12p . Without loss of generality, we assume / to be a real valued function. The first 
identity (I3.25P follows readily by noticing that 



|Vs2(/®eJ(^,0)P = — 



27r 



m 



sin' 9 



urn' + If 



The proof of ( ]3.26p requires more calculations and application of integration by parts sev- 
eral times to take care of some cross products appearing in the integral form of the norm. 
Without loss of generality, for a fixed m e Z, we can assume / G span {Q™ : n > \m\} 
because this subspace is dense in W^. Observe that 



/(O) = f (0) = fin) = fin) = 0. 



Since 
1 



m 

+ 



2 /-I- M \ 



(A.l) 



|/'(0)|2sin^d^ 



16 



/2(e) sin^de. 



it is sufficient to analyze the term containing the Laplace-Beltrami operator: 



|A§2(/ (g) e„)||^o 



^0 



|A§2(/®e^)(e,0)|^sineded0 



m 



. 2./W + -^l(sme/'(e))' 
sm Q sm ' ^ 



sin^de 



de 



(sine/(e))'|^de 



sin Q Jq sin 6 

,/ de 



-2m' / f{9) ( sin Of (9)) =: m'h + h- 2m' I,. (A.2) 

' sm 9 



Using (lA.ip and integration by parts 

"cos^ 9 



sin 9 

" cos2 9 . ^, 



\n9)\' +M\f"{9)\' + cos9{\f\9)\')' 



d9 



sin 6* 



\f\9)\'d9+ / sin9\f"{9)Yd9 



d9 



\f\9)\' sin9 A9 



\f\e)V—n+ / ^^r^e\n9)\'d9. 



(A.3) 
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Proceeding similarly, we derive 

Jo \sm^9 J 

Jo sm6' Jo sm^fc' smt^^ 

= - r\m\'^,- r rm^de. (A.4) 

Jo sin6' Jo sine* Jo sin.-" 6 

Inserting (lA.3p - (lA.4p in (lA.2p . we obtain the desired result. □ 



Proof. (Corollary [321) For |m| > 1, fICTll follows from ([325]). 



The inequalities 



1 

16 ^ 4 



9m^ 1 

+ (m^ - 4m2) < < 3m^ ^ + (1 + 2^^) < V|m| > 2 



with (I3.26P imply the first inequality of (13.300 . For |m| > 3 the second inequality of (13.300 
is simply a consequence of the inequalities 

4 2 

- 4m^ > m'^/2 > — , 1 + 2m^ > > — . 

^ 6 ' - 6 

The case |m| = 2, has to be analyzed separately since one of the crucial terms, the third 
term in (I3.26p . vanishes: Using (I3.26p . we obtain 



2 > 9 



d^ 



sin 9 



9 



AO 
sin 9 



\f\e)\\inede. 



(A.5) 



As before it suffices to consider / to be real valued and that / G span {Q^ : n>2\. Note 
that 



sm U In. sm 6 



dd. 



(A.6) 



nef^de 



de. 



Applying integration by parts to the second term and using flA.lO we obtain 

cos 9 

sin^ 6 

Notice that for 9 E (0,7r/2], 

sine|logtan(^/2)| < 2tan(^/2)| logtan(e/2)| < 2e"^ < 1. 

By symmetry, we can extend this bound for any 6 G (0,7r). With the help of (lA.Sp and 
the inequality 2ab < a'^ + b'^, from (IA.7P we obtain 



(A.7) 



(A.8) 



nef-^do < 



sin"^ 9 







1 

< - 
- 2 



df) 

f{0)f\e) ^ + 

smt^ 



d^ 



f{e) — + 2 \f 



, d^ 
sin 9 



1 

+ 2 



sm 
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Inserting (]A.9|) in ( 1A.6|) we easily derive 



'sin^6' 2 Jq ' ^"'sinO ' Jq " sin^ ' 2 Jq ' ^''^sin''6' 
and therefore 

< 3 Ti/wp^- (A.io) 

From ( lA.SP and (lA.lOp . we obtain 

PTT JQ l-TT 1/] f-n 

611/11?., > 54/^ m^g^^l \fv^\'--^^^ 

> i6(3rmii^+2rimp^) 

V Jo sme Jq smej 

+4 r\f'm'^+ r\n9)\\m9d9 
Jo sm 6* Jq 

PIT J/1 /-TT 1/1 pn 

> 16/ |/(^^)p-^ + 4/ \rm^—+ \ne)\hineAe. 

Jo sm (7 Jo smt7 Jq 

Hence the inequahties in f l3.30p hold. □ 

Proof. (Proposition 13.41 ) Denote by F the maximum circle in S^, parametrized by 

g(0) := (sine,0,cose). (A.ll) 

Given /° : F — )■ C we denote / = /° o q : M — )• C. The norm in the Sobolev space ^^(F) 
can be then defined with the help of q and (13.321) : 



\\r\\H-{r) ■■= \\J \\h: 



# 



The second ingredient we will use in this proof is the trace operator 7r which can be 
shown to be continuous from '^'""'"^''^(S^) onto H^{T) for all r > (see [71 |20j for a proof 
of this result in M"; the proof can be easily extended by using local charts of the unit 
sphere and the equivalent definitions of the Sobolev spaces involved). 

Given / G W^, consider the mapping 

V^f := v^(7rF°) o q, F° := (/ ® e„) o p-\ 

Observe that F° G 'H'"(S^) and that actually / = Vmf, that is Vm is simply the identity 
operator. Moreover, 



"Pm/ll/fJ < V27r||7r||^r+l/2(§2)^jyr-(r)||-F°||^r+l/2(g2) - V 27r 1 1 7r 1 1 ^r + 1/2 (g2 )^ jyr (p) 1 1 / 1 1 ^;(^r^+l/2 , 
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where i|7r||-H'-+i/2{§2)_5.jyr(r) is the continuity constant of 7r as a hnear operator from 
^r+i/2^g2^ onto H''(T). Hence we obtain (Km . 

Since = L^j^, f l3.35p is clear whereas fl3.36p for m = follows directly from f l3.25p and 
(1333]). 

Finally, if m 7^ using /(O) = (see Remark I3.3p . we observe that 



d^ 

\m\'- 



sin 9 



Tv/2 



sin 9 



d^ < 



-/^ V9 
sin 9 



i/(orde 



d^ 



Proceeding similarly, but using now that /(vr) = 0, we can bound the integral in (7r/2, vr) 
and hence conclude that 

r\9 

\fm'—-n<c i/(orde 



sin 9 



Equation ( ]3.29p now yields that 



V5 



lyi < VII/IU,^ <C(l + |m|)||/||^i. 



□ 



Proof. (Proposition 13.51 ) For the sake of a simpler notation, we will assume without loss 
of generality that / is a real valued function. We consider the functions 

f 1, if^<7r/6, 

MO) ■= i 2--, if ^^G [7r/6,7r/3], 
vr 

0, otherwise. 

Obviously, < ipj < 1 and actually {ipi,ip2,ip3} C H^{0,tt) form a partition of unity. 
Thus 

3 3 

r fi9) d9 = J2 r fiO)MO) =: (A.12) 
Jo „_i Jo „_i 



Clearly 



' fi9)ip2i9)d9<2 I ' f{9)sin9d9. (A.13) 

7r/6 J 7r/6 



On the other hand 



h = f\9)^,{9)d9 = f{9)^,{9)9'^^'' - r'\f{9)^m)'0d9. 
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The first term is zero, since <y9i(7r/3) = 0. For the other term we first note that 

6 



Thus, 



7r/3 



||V^'l|U°°(0,7r) < 

27r 



^ 1.209, V^e[0,7r/3]. 



< 



< 



< 



3^3 
4 

7! 

4 

7! 



2 r^'\fie)f{9)\sm9de+ fie)Me)\smede 

.Jo Jo 

7r/3 7r/3 

2/ \fi9)f{9)\smede+ [ f{9)sm9d9 
Jo Jo 

7t/3 \ 1/2 / /.7r/3 \ 1/2 



sin 9 3^3 



(2/(^)/'(e) + f{9)cp[{9))-^ sm9d9 

Sill u 

7r/3 



f\9) sm9d9 





f7r/3 



sin^d^ 







+ / f{9) sm9d9 
Jo 



(A.14) 



Proceeding in a similar way, we can prove 



h < 



f{9) sm9d9 



27r/3 



1/2 



\f{9)\ sm9d9 



27r/3 



1/2 



sin^d^ 



27r/3 



Substituting ( ]A.13l) -( p^.l5l) in ( ]A.12I) . we obtain the desired result. 



(A.15) 
□ 



Recall that if Q is open and F C f2 is a Lipschitz curve, then [HI Theorem 1.6.6] 



l7rM||i2(r) < Co||M||L2(f7)||M||Hi{Q)- 



It is easy to see that Proposition 13.51 can be alternatively derived from this result by taking 
i7 = S^, F the circle parameterized by flA.lip and u = f o e^- 



B Proofs of Lemma 14.51 and Proposition 14.6 



In this section we prove several auxiliary results, most concerned with convergence esti- 
mates of the interpolants q™,gi' introduced in fl4.7p - fl4.9p . in the spaces W^, which are, 
up to our knowledge, new or sharper than that found in the literature (see for instance 
[3] or [Til §6.6] and references therein). 

The first result is an inverse inequality for and which holds in VF^, uniformly in 
m and N . 
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Lemma B.l There exists C > such that for any m & qn & ©at U with 5'Ar(0) = 
g^(iT) = zfm^O, 

\\9n\\w^ < C{1 + |m|)A^||^jv|Uo . 

Proof. We first prove the result for even m. Let qn G and take m' = 2 if 5'Ar(0) = 
QNij^) = and m' = otherwise. Note that 

fifiv G span {Q™' : n = m', . • • , A^}, 

(this property is not vahd for m' > 2). Since {Q^'}^=rn' orthonormal basis of 

^sin = W^m. we obtain 



/ L^o 



N 



Qn= 



n=m' 



Hence 



I l|2 

\9n\\wi^ 



N 



N 



n=m' 



This concludes the proof for |m| = 0, 2. For even \m\ > 2, we apply fl3.25p of Proposition 
IS.ll to derive 

I m I I m I I m I 

\\9n\\w^ < —\\9N\\wi<—iN + ^)\\gN\\q.^ = —{N+l)\\gN\\wg- 

The proof for odd m is similar: Since = span {Ql^ : n = 1, . . . , N}, as above the result 
can be derived first for m = 1 and using again f l3.25p of Proposition I3.H the result follows 
for all odd values of m. □ 



Let 



' n 1 



fmin)= I fie)Q^ie)smede, 

10 



(B.l) 



\m\<n<N-l 

with the convention that, for |m| > N, s^f = 0. Clearly s^ is just the orthogonal projec- 
tion onto span {Q™ : n = \m\, . . . , iV — 1} in for all s. Moreover sj^f G 
s^+V e a^iid therefore 



o ^21+1 _ 2e+l 



(B.2) 



Lemma B.2 For all r > s and N > 2, 



Moreover, if m 0, 



ii/-s^/iu. <(iv+ir'^ 



s^/(0) = s^/(7r) = 0, 



(B.3) 
(B.4) 
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and 



for all r > 1. 



max 



{|sU(0)-/(0)|,|sUW-/W|}< 



l-r 



(2r- 2)1/2 



(B.5) 



Proof. Take A^' = max{A^, |m|}. Then 



2r,V / x|2 

n r 



n>N' 



n>N' 



l\2s-2r.: „||2 



which proves the first result. Using fl3.3p -f l3^ . 

Qn (0) = Q'^in) = 0, for m ^ 0, 

and hence (1B.4P follows. 

For m = we note first that 

Qn(0) = (r^+i)'/' = (-l)"g°(vr). 

Then 

14/(0) -/(0)| < $^|/oW|g°(0) = ^(n + i)^/^|/o(n)| 



n>N 



< 



< 



E("+5) 

n>Af 



n-(2r-l) 



.1/2 



.1/2 



n>Af 



(2r- 2)1/2 
where, for a > 1, we have apphed 

j:{n+ir< / (x+i)-"dx<-i^iv-". 

Proceeding in a similar way, we can bound |s^/(7r) — /(vr)| and hence we obtain (IB.Sp . □ 



^ the 



Our next result gives us information about the distribution of the nodes {6j}j^Q, 
nodes of Qff in the elevation angle. Roughly speaking, this distribution turns out to be 
quasi-uniform. 



Lemma B.3 For N > define 

. 0-l/4)7r 



O'N/2 



N+1/2 
TT _ 57r/8 



2 
TT 



Af+1/2 
0-N-j 



'^N/2 



Ar+l/2' 
TT I 57r/8 
2 Af+1/2' 

7T — b 



N- 



for even N , 

j = m,...,iv-i. 
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Let {Ojjf^Q be as m (EH]). Th 



en 



Proof. Note that for even N, 6^/2 = arccos^7v/2 = arccosO = 7r/2 and therefore there is 
nothing to prove in this case. On the other hand, 

arccosCAT-j+i < Oj = arccos?7Ar_j < arccosCAf-j 

where {Cjj^i are the zeros of P/v, the Legendre polynomial of degree N (see (13. 3p ). 
We can apply well-known results on the asymptotic distribution of these roots (see for 
instance [III (3.1.2)-(3.1.3)]) to deduce that for 1 < j < N/2, 

(j - l/4)7r ^ ^ , (j + l)7r , 

^ N+l/2 ^ a^ccosCiv-i+i < Oj < arccosCiv-i < ^ ^ = bj. 

(See also [271 Theorem 6.21.3] or [3l Theorem 2.4].) For j > N/2 the result follows from 
the first part of the proof and the symmetric disposition of the nodes. □ 



Let us introduce the Gauss-Lobatto quadrature rule 

N 1 

which is exact for any polynomial of degree up to 2 — 1 [HI [H] . We point out that the 
weights of this formula satisfy (see for instance [HI Section 3.5], [3^j Corollary 2.8]) 

2 



N{N + 1] 

and 



0<a;,<^sin^„ i = l,...,iV-l (B.6) 



with C independent of A^. 
Consequently, the quadrature rule 



CnI ■.= ^^jf{ej)^ f{9)sm9d9= /(arccoss)ds 
j=o Jo J-i 

is exact for any G in'2Ar-i (recall that 9j = arccos?7Ar_j). Hence, 

Cn{\pn\^) = WpnWI^ , e U D^_2. 

sm 

If Pat G U (see (12. 6p ). the rule fails to be exact for Ipat^ G ^2N- However, we 

have instead the bound 

/ \pM\'^^^dd9<Cj2^n\PN{0n)\^ PNeB>%UB>%_, (B.7) 

n=0 
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with C independent of N and pn- This bound, which is a consequence of a well known 
stability result for Gauss-Lobatto formulas (see [HI Theorem 3.8.2]), will play an essential 
role in what follows. 

We are ready to prove the first inequality in Lemma 1431 In fact below we prove a slightly 
more general result (see Corollary 13.21 which gives an equivalent norm for W^). 



Lemma B.4 There exists C > such that for all N > 2 and for any f E C°[0, vr] with 
IKgi/ll^L ^ ^fll/II^L + ^ll/'II^L + 4(1/(0)1 + 



with C is independent of f and N . 



Proof. Without loss of generality we assume that / is a real valued function. Recall that 
from the definition of q7v„i, 



q#,gi/(^.) = <tNj{o,) = f{e^), j = o,...,iv. 



qS/(^.) = q^,gi/(^.) = /(^.), J = l,...,iV-l, q^^+,V(0) = q^i /(vr) = 0. 



21+1 . 



,21+1 



Then, using (IB.7P we obtain 



N 



N 



n=0 



n=0 



where Cm is independent of / and A^. Note that Cm = 1 for m odd (that is, for q^gi/) 
or for m even with q^gi/ € ID'^„]^. From Lemma [B. 3 1 we have 

9ne[an,bn], with ^ < 4{N^+l/2) ^ ~ ""^ - 1^' n = l,...,N-l. (B.8) 



Hence from ( ]B.6p and ( IB. 81) we conclude that 

2Crn 



[q^,gi/]^(^)sin^d^ < 



< 



N{N + 1) 

2Cm 

iv(ivTi) 



(f(0) + f(vr))+C 



N-l 



^^/2(^„)sin^„ 



n=l 



c 



N-l 



(f (0) + fin)) + - J2ibn - an)/'(^n) siu 



n=l 



(B.9) 



We point out that the above sum can be understood as a composite rectangular rule 
applied to f'^{6) sinO. Using the well-known bound (with c G [a, b]), 



g{e) Ae-g{c){b-a) 



<ib-a) \g 



d9. 
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and noticing that 



J2 / \9{0)\de<2 / \gmdO, WgeL\0, 

n=l 



we deduce 

N-l 

- an)f{9n) sin 9r. 



n=l 



N-l 



n=l 
N-l 



^-1 fbr. 

n=l 



<^(6„-a„)/ \{f{e)smey\de + 2 f{e)smede 



An 

< — 

C 

< — 

- N 



2 / f (9) sin Ode 
Jo 



1 f{e)\cose\de+ / 2\f{e)f'{e)\sinede 

Jo 

/•TT / /»7r \ 1/2 / /»7r 

' /2(0)d^+( / f{e)sinede\ ( / |/'(0)psin0d0 

L"'o \Jo / \Jo 

/•TT 

+2 / f{e)sin9de 
Jo 

pTT /-ill / pTT \ 1/2 / /.TT 

<C"y /2(^)sin^d^ + -^( / /2(0)sinede^ ' ' 



1/2- 



1/2 



I If {6)1^ sine de 

lo / \Jo / 

where we have apphed again (IB.SP and used in the last step Proposition I3.5[ Finally, the 
inequality 

2ab < Na^ + N'^h^ 

proves that 



N-l 



J2ibn-an)fie^)sine^^<C 



n=l 



f{e) sin^d^ 



/o 



|/'(^)|2sin^d^ 



Substituting (IB.lOp in f lB.9|) . the desired result follows. 



(B.IO) 

□ 



Remark B.5 Inequalities such as those proved in Lemma IB.4I have been used in the 
literature for studying the convergence of polynomial interpolants in weighted Sobolev 
spaces (see for instance [TTj or [3] and references therein). 

Note that with the change of variables s = cos 6', Lemma IT5] implies that for / G 1, 1] 
with /(±1) = 0, and if pn is the polynomial of degree N which interpolates / at Gauss- 
Lobatto nodes {?7j}^0' have the bound 



\pN{s)\'ds<Ci / \f{s)\'ds + 



\nsni-s')ds 
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with C independent of N and /. This bound is an improvement of Theorem 4.1 in [3] 
where a similar result is proven but without the weight function (1 — s^) in the second 
integral. □ 



Next we prove Proposition 14.61 for the s = case. 
Proposition B.6 For r > 1 there exists > so that 

iiq^,gi/-/iko <aiV"ni/ik^- 

Moreover, the estimate above holds also for r = 1 if m ^ 0. 



Proof. Assume first that m ^ 0. We recall from Remark 13.31 that any / G 

is continuous and vanishes at {0,7r}. From (1B.2I) . Lemma IB. 41 (see also Lemma 14. Sp . 

Proposition 13.11 and ( IB. 31) of Lemma [B. 21 we deduce that for r > 1, 

gi/-/lko = ||q™^^(/-s-/)-(/-s-/)|Uo 

< C[\\f -s^fWwo^ + N-'Wf -s^fWw^] <CN-^\\f\\w^. 



If m = 0, we have instead 



qkgi/-/ll< < C(\\f-s%f\\^o + N~' ||/_sO^^||^^^ + |(^_gO^^)(0)| + |(/-sU)WI 



Using ( IB. 31) as before, and ( IB.SP for taking care of the pointwise errors at {0, tt} (and here 
using the restriction r > 1) we conclude that 



HN,,J-f\\wS<CN-^\\f\\ws- 



□ 



To derive convergence estimates in W.^ we require additional technical results. To this 
end, it is convenient to introduce the norms 



sin^ 



nl/2 



\f{e)\'sm^ede 



and the associated spaces L^. ^.. Below we will use these norms for k = —3, —1 and 2 



Lemma B.7 For all N >2 and for all eB%U lED^^i, 



L2(0,7r) 



< 



sm 



1/2, 
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Proof. In this proof we will apply the fact that the quadrature rule 

N+l ^ 



is exact for any / G '^2N+2- Notice that, by hypothesis, \pn\'^ ^ "^tN- Hence, 

N+l 



N+l 



< 



sin (2]^) 



vr 



2„ „2 
L2 



fB.lll 



We point out that, in the usual notation of the theory of interpolation of Hilbert spaces, 
we have 

4^in=[^^(0,vr),LL.],/, 
(see for instance [281 Lemma 23.1]). Hence, with (IB. lip , we obtain the desired result: 



Pjv|li2(o,^) < [sin(2#^)]"'|bjv|li2 , Vp^veo^^uro^.i. 



□ 



Lemma B.8 There exists C > so that for all f G with \m\ > 2, 



^ sin ^ 



^2 . + 



1 

N 



-3 + 



L2 

sin 



J 



(B.12) 



Proof First, consider the case of even m so that q^v^gi/ = 'iN,gif- Then, it is easy to 
check that 

{q%,,jm/sme = q%Jf/sm)ie). 

(Note that in this case /(O) = /(tt) = 0.) We can proceed as in the proof of Lemma [3.41 
see flROjl -f lRTOj) . and derive 



7V-1 



llAf,gl/llL2 



n=l 



= ||q^,gi(//sin)||i,_^ <C5^(6„-a„)(//sin)'(0„)sm^„ 

< c 

< 2C' 



IJO 



^'(^^i^+i^ / |(//sin)>)|%in^^de 



L ^0 



sm N'^ Jq sm 9 Jq sm 
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This proves the result for the even m case. 

To prove the result for the odd m case, observe that now q^gi// sin G ^%^2- However, 
if we try to apply the same ideas as before, unlike the previous case, we have to take 
care also of the pointwise values at and vr. Indeed, following again ( ]B.9[) -( lB.10[) . and 
noticing that Cm = 1 in ( ]B.9p . we derive 

N-l 

+C^(6„ - a„)(//sin) (6'„)sin6'„ 



< 



n=l 

2 



N{N 
+2C' 



-{\{q%,,J/ sin ) (0) r + I (q°Ar,gi// sin ) (tt) 



f\0) 



de 



+ 



IJO 



sine 



+ 



sin^ e iV2 Jo 



de 

sin^ 



(B.13) 



It remains to bound 

RnU) := iv(iV+ 1) (l(q°^,gi//^i^)(0)l' + l(q°^,gi//sm)(7r)r) 
4 



< 



A^(A^ + 1) "q^.gi'^/^^^ IIl-CCtt)- 
We proceed as follows: Assume first that / is even, that is 

f\7,-e) = f{e). 

Then, so is (f^^^J so that 

L(7V-2)/2j 

(qAr,gi//sin)(6') = ^ a2nCos2ne 

n=0 

for suitable coefficients a2n- Hence, by applying the Cauchy-Schwarz inequality we obtain 



|q^,si// sin II 



L°°(0,7r) 



< 



< 



[{N~2)/2\ 



n=0 



\0i2r, 



IT 



[(N~2)/2\ 



TT 



i"oi'+ Yl 



\a2n\ 



n=l 



2J 



X 



N N. 



TT 



q^ i// sin 11^2 



(O.tt)- 



Recall that (1% ^if / sin G 



_i. Then we can apply Lemma [B. 71 to get 



N 



|q^,gi//sin|Loo(o,.) < ^gin(^/((2iV + 2))"'^^'S^^/'^''"^?in' 
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and therefore 

RnU) < 



7r{N + 1) sin(7r 



/P^l|q°..,,//sin||i.^ =: c(iV) ||q^,,/||i.^_^. 



We point out that 



< c(N) < c(2) = — < 1, ViV > 2. 



Using this bound in ( 1B.13|) . we obtain 
2C 



glJ II L2 



< 



1 - c(2) 



, d^ 
sin 6 



(B.14) 

For odd f, that is for functions with f{6) = —/{it — 6), we can proceed similarly and show 
that a similar bound to flB.14p also holds in this case. 



We have proven the result for the subspace of even and odd functions. To extend the 
result for arbitrary / we note first that 

f = fe + /o, 

where fe{0) '■= |(/(^) + /(^ — ^r)) and fo{0) = ^{f{0) — f{6 — tt)) are even and odd parts 
of /. We deduce that 



I^Af.gl/llia 



I^Af.gl/ellia + ||qAr,gl/o|lL2 



< c 



c|Il2 , + II /o 11^2 



ll/<"^ 
(ll/< 



jy2 VII^c|Il2 



c 



3 + ii/oiii2 3 + ii(/eyiii2 ^ + ii(/oyiii2 j 

sin sin -'- sin -'- 

1 



+ 



2 

L2 



+ \\f\\h J 



Thus the result (lRl2l) follows. 

We are now ready to prove the second inequality in Lemma IB.4[ 



□ 



Corollary B.9 For all \m\ > 2, 



71 H 



Proof. Note that 



|q^,gi/ll 



< m 



/||2 



"'q^,gi/lli2 +2||(q^,gi(/-s^/))'||^, +2||(s^/) 

t^in J- sin 



III 



=: /1 + /2 + /3. 
We bound the first term using Lemma [B. 8 1 



Ji < Cm^ 



f'l|2 



< c 



Zl + 



Ar2 
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For the second term, using first Proposition 13.11 and combining the inverse inequahties 
stated in Lemma IB. II with Lemma IB. 21 and Proposition IB.6t we obtain 

h < 2||q^,gi(/-s-/)||J,, < CiVlq-g,(/-s-/)||J,, < C'||/f^. < 
where in the last step we have used 03.291) of Corollary 13.21 We bound the last term as 

/||2 ^ oil /Il2 ^ ^ II ^l|2 



^3 < 2||s^/||^. <2||/||^^. <-,„„^. 
Thus we obtain the desired result. □ 

Finally, we prove Proposition 14.61 for the s = 1 case. 
Proposition B.IO For all r >2 there exists Cr > so that 

\hNsf-f\\w^<CrN'-'-\\f\\w^. 

Proof. Note that 

\hN,,lf - fWw^ = hN,M-^Nf)-if-^Nf)\\wl 

< hNjf-^Nf)\\wi^+N'-'\\f\\w^ 

where we have applied flB.2p and flB.3p of Lemma IB. 21 

It remains to bound the first term. For |m| < 1, we make use of Lemma IB. II and 
Proposition IB. 61 to obtain, 

l|q^,gi(/-s^/)lk^ < C'iV||q-g,(/-s^/)||^o <CAr[||q-^J-/i|^o +||/-s^/|U^^ 

< C'N'-^WfWw^. 

For |m| > 2, this approach is not valid since the inverse inequalities of Lemma fB . 1 1 contains 
m as a penalizing term. Thus, we follow a different approach and use Corollary 13.21 and 
IB. 91 to obtain 



< C^(ll/-s^/|k^ + V^iV-i||/-s^/||H..) <C"iVi-ni/|k;;. 
Hence the desired result follows. □ 
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